### 5.1 Editando alguns parmetros grficos
e.saligna =  read.csv("esaligna.csv", sep = ",", header = T)
head(e.saligna)
str(e.saligna)
e.saligna$talhao = as.factor(e.saligna$talhao)
par(cex.axis=1.5)
par(cex.lab=1.5)
plot(e.saligna$dap,e.saligna$ht,  bty="l", tcl=0.2, cex = 1.5, xlab = "Dimetros  altura do peito (cm)" ,  ylab = "Altura (m)")
mtext("rvores de Eucalyptus saligna", line=0.4, col="red",cex=1.5)
### 5.2 Dois grficos juntos
par(mfrow=c(1,2))
boxplot(e.saligna$dap~e.saligna$talhao, ylab= "Dap", xlab = "Talhao", col= "grey")
text(x=0.75, y = 22, "a")
media = aggregate(e.saligna$dap,by=list(e.saligna$talhao), FUN=mean)
class(media)
desvio = aggregate(e.saligna$dap,by=list(e.saligna$talhao), FUN=sd)
g.talhao = data.frame(media,desvio)
g.talhao = g.talhao[-3]
colnames(g.talhao)[1] = "talhao"
colnames(g.talhao)[2] = "media"
colnames(g.talhao)[3] = "desvio"
head(g.talhao)
str(g.talhao)
plot(g.talhao$media, xaxt="n", ylim=c(0,25), xlab="Talhao", ylab = "Dap")
axis(1, at=c(1:6),label=c("16","17", "18", "22", "23", "32"))
segments(c(1:6),g.talhao$media-g.talhao$desvio, c(1:6), g.talhao$media+g.talhao$desvio)
text(locator(1),"b")

### Exercicio 5.3
ex = read.table("exercicio3.csv", sep = ",", header = T)
head(ex)
X11()
par(mfrow=c(1,2))
par(cex.axis=1.0)
par(cex.lab=1.2)
plot(ex$x1,ex$y1, xlim = c(0.45,2.2),ylim=c(0,3.0),tcl=0.2, cex = 1.0, bty = "l", pch = 17, ylab = "Euclidian distances", xlab= "Log(Patch size)(ha)", )
segments (0.44,2.18,2.2,1.5, lwd = 2.0)
text(locator(1) ,"a", cex=1.5)
par(tcl=0.5, bty="l")
boxplot(ex$x2~ex$y2, ylim=c(-0.1,3.1), xaxt = "n", yaxt="n", tcl = 0.2, cex =1.0, outline=F, bty="l")
par(cex=1.0)
axis(2,at=seq(0,3.2,0.5),tcl = 0.2) 
axis(1,at=c(1:6), tcl = 0.2,label=c("Small","Medium\nEdge","Medium\nInterior","Large\nEdge","Large\nInferior","Control"))
text(x=c(1:5), y=rep(3,5), c("*", "*", "**", "*", "***"))
text(x=6.5,y=3.1, "b",cex=1.5)